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We calculate the equation of state of isospin-symmetric nuclear matter in the three- 

■ loop approximation of chiral perturbation theory. The contributions to the energy per 
CN ! particle E(kf) from one- and two-pion exchange diagrams are ordered in powers of the Fermi 

momentum fcj (modulo functions of kf /m n ). It is demonstrated that, already at order O(kj), 

■ two-pion exchange produces realistic nuclear binding. The underlying saturation mechanism 
is surprisingly simple (in the chiral limit), namely the combination of an attractive fc^-term 
and a repulsive fc|-term. The empirical saturation point and the nuclear compressibility 

\ K ~ 250 MeV are well reproduced at order 0(kf) with a momentum cut-off of A ~ 0.65 GeV 



f> 

which parametrizes short-range dynamics. No further short-distance terms are required in 
our calculation of nuclear matter. In the same framework we calculate the density-dependent 
asymmetry energy and find Aq ~ 34 MeV at the saturation point, in good agreement with 
the empirical value. The pure neutron matter equation of state is also in fair qualitative 
agreement with sophisticated many-body calculations and a resummation result of effective 
^ ! field theory, but only for low neutron densities p n < 0.25 fm~ 3 . 



PACS: 12.38.Bx, 21.65.+f 
Keywords: Effective field theory at finite density; Nuclear matter equation of state; Asymmetry 
energy, Neutron matter. 



1 Introduction 

One of the basic problems in nuclear physics has traditionally been to develop a microscopic 
understanding of the nuclear matter equation of state and the properties of finite nuclei in terms 
of the underlying free nucleon-nucleon interaction. The present status is that a quantitatively 
successful description of nuclear matter can be achieved, using advanced many-body techniques, 
in a non-relativistic framework when invoking an adjustable three-body force 0. Alternative 
relativistic approaches treat nucleons as Dirac-particles. The simplest one in this category is 
the <7u;-mean field model of Serot and Walecka 0. It describes nucleons as Dirac-quasiparticles 
moving in self-consistently generated scalar and vector mean fields. With two couplings adjusted 
to the empirical saturation point, the nuclear compressibility comes out however much too large 
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(by a factor of 2). Refinements of relativistic mean field models which include additional non- 
linear terms with adjustable parameters, are nowadays widely used for the calculation of nuclear 
matter properties and finite nuclei |3|]. The more basic Dirac-Briickner approach of Brockmann 
and Machleidt solves a relativistically improved Bethe-Goldstone equation with a realistic 
nucleon-nucleon interaction parametrized in terms of one-boson exchange potentials. 

During the last decade, a novel approach to the NN-interaction based on effective field theory 
(in particular chiral perturbation theory) has emerged [^, |], [7|, §]. The key element is a power 
counting scheme which relies on the separation of long and short distance dynamics. The meth- 
ods of effective field theory have recently been applied to systems at finite density in refs. |9|, [T(| . 
These works deal mainly with the effective range expansion (and generalizations thereof) of the 
equation of state of a dilute many-fermion system. The complete resummation of in-medium 
multi-loop diagrams for a system with an unnaturally large scattering length (such as neutron 
matter) has been achieved in ref. (in the limit of infinite space-time dimensions). Lutz et al. 
[ |l2|j have calculated the equation of state of nuclear matter from chiral pion-nucleon dynamics. 
The contributions of one- and two-pion exchange diagrams to the energy per particle have been 
evaluated up to fourth order in small momenta. An essential ingredient in their work is an addi- 
tionally introduced attractive zero-range NN-contact interaction. The saturation point obtained 
in this calculation depends very sensitively on the numerical value of g\/2 + go + g\, the strength 
of this NN-contact interaction. The dynamical origin of the important attractive NN-contact 
interaction remains unexplained. Furthermore, the two components of the contact interaction, 
the one proportional to g\/2 and the other one parametrized by a coupling <?q + #i, have not been 
treated on equal footing in ref.||12|| when going to higher orders in perturbation theory. 

In its starting point our approach to the nuclear matter problem is closely related to the work 
of ref. [T^j . We also use the three- loop approximation of chiral perturbation theory at finite density 
to calculate the contributions of one- and two-pion exchange diagrams to the energy per particle 
of isospin symmetric nuclear matter. In a systematic non-relativistic expansion we compute all 
terms up-to-and-including fifth order in small momenta (these are the Fermi momentum kf and 
the pion mass m^) and thus go one order beyond the work of ref.[|12|. A major difference in 



comparison with ref. |12j is that we do not introduce any zero-range NN-contact interactions with 
adjustable strengths in addition to the chiral pion-exchange terms. Instead we use a momentum- 
space cut-off A to regularize the few divergent parts associated with the chiral two-pion exchange, 
and we keep the power divergences which are specific for cut-off regularization. Their contribution 
to the energy per particle is equivalent to that of a zero-range contact interaction parametrized 
through the scale A, which then represents the short- distance dynamics. This relationship is very 
useful in order to understand the saturation mechanism underlying the chiral two-pion exchange, 
as we shall discuss. A further feature of our treatment is that the properties of isospin symmetric 
nuclear matter (saturation point and compressibility) and of isospin asymmetric nuclear matter 
(asymmetry energy and equation of state of pure neutron matter) are given in terms of only one 
adjustable parameter, namely the cut-off A. It is then highly non-trivial to satisfy all known 
empirical constraints by fine-tuning just this one single scale A to a physically sensible value. We 
demonstrate that this is indeed possible. The present chiral perturbation theory calculation allows 
to systematically explore the role that chiral symmetry plays in the nuclear matter problem. 

Our paper is organized as follows. In section 2 we perform the chiral (or small momentum) 
expansion of the nuclear matter equation of state up-to-and-including terms of order 0{kj). We 
present analytical results for the energy per particle E{kf) as given by l7r-exchange, iterated 
l7r-exchange and irreducible 27r-exchange diagrams. Then we discuss the saturation mechanism 
and the properties of our nuclear matter equation of state. In section 3 we calculate the density 
dependent asymmetry energy A(kf) in the same framework and discuss the corresponding results. 



Section 4 is devoted to the equation of state of pure neutron matter as it follows from the same 
set of pion-exchange diagrams. Finally, section 5 ends with a summary and an outlook. 




Fig. 1: Equation of state of is o spin symmetric nuclear matter. The full line gives the parametriza- 
tion eq.(l) with a = 5.27 and (3 = 12.22. The dashed line stems from the many-body calculation 
ofref.ftrg. 

Before presenting results of our (perturbative) diagrammatic calculation we like to draw at- 
tention to the following fact. A simple but nevertheless realistic parametrization of the energy 
per particle of isospin symmetric nuclear matter is given by the following expression, 

op p p 

with kf the Fermi momentum related to the nucleon density in the usual way, p = 2kj/3ir 2 , and 
M = 939 MeV the (free) nucleon mass. The first term in eq.(l) is the well-known kinetic energy 
contribution of a Fermi gas and the powers of M in the other terms have been chosen such that a 
and /3 are dimensionless. The two-parameter form eq.(l) has generically a saturation minimum if 
a, [3 > 0. The interesting feature of this parametrization is that once a and f3 are adjusted to the 
empirical nuclear matter saturation point, fc/o = (1.35 ±0.05) fm _1 and Eq = E(kfo), the nuclear 
matter compressibility ]13|, |i~4"f , 



k = * d2E ^ 



(250 ± 25) MeV , (2) 



comes out correctly. For example, adjusting to the saturation point obtained in the sophis- 
ticated many-body calculation of the Urbana group [fl5|l , po = 2/cj /37r 2 = 0.159 fm -3 and 
Eq = —16.0 MeV, one finds a = 5.27 and (3 = 12.22, and the compressibility is predicted as 
K = 236 MeV. In Fig. 1 we compare the density dependence of E(kf), given by eq.(l) (full line) 
with the equation of state resulting from the many-body calculation of ref.|T5|] (dashed line). One 



observes that the relative deviations of both curves do not exceed 10% up to quite high densities, 
p = 1.0 fm -3 . Keeping in mind that eq.(l) provides a realistic parametrization of the nuclear 
matter equation of state will be very useful in order to understand the saturation mechanism 
underlying the chiral two-pion exchange. In the chiral limit (m n = 0) and truncating to order 
O(kj), the equation of state following from one- and two-pion exchange will be precisely of the 
form, eq.(l). 

2 Chiral expansion of the nuclear matter equation of state 

The tool in order to systematically investigate the consequences of spontaneous and explicit chiral 
symmetry breaking in QCD is chiral perturbation theory. Observables are calculated with the 
help of an effective field theory formulated in terms of the Goldstone bosons (pions) and the low- 
lying baryons (nucleons). The diagrammatic expansion in the number of loops has a one-to-one 
correspondence to a systematic expansion of observables in small external momenta and meson 
(or quark) masses. For the problem under consideration, the external momentum is the Fermi 
momentum kf, related to the nucleon density p = 2ky i /3ir 2 . The relevant observable is the energy 
per particle E(kf), i.e. the ratio of energy density and particle density, with the free nucleon 
mass M subtracted. Consequently, the equation of state of nuclear matter as given by chiral 
perturbation theory will be of the form of an expansion in powers of the Fermi momentum kf. 
The expansion coefficients are however non-trivial functions of kf/m^, the dimensionless ratio 
of the two relevant small scales inherent to the problem. Note that at the empirical saturation 
point, kfo ~ 2171^, so that both scales are of comparable magnitude for the densities of interest. 
This also implies that pions have to be kept as explicit degrees of freedom. At the saturation 
density pionic effects cannot be properly accounted for just by coefficients of local NN-contact 
interactions which form the basis of effective range expansions at finite density ||10| . 

Naive chiral power counting (which is basically a counting in terms of mass dimension) suggests 
that closed vacuum diagrams with L loops (which represent the ground state energy density 
in diagrammatic language) give rise to a contribution to the energy per particle of the form 
E(kf) = k"j L ~ l J r L(kf/m 7T ). However, the two-nucleon system is known to provide exceptions to 
the naive counting rules by the so-called iterated one-pion exchange. In that case the relevant 
energy denominator is a difference of nucleon kinetic energies, and this implies that the large scale 
factor M (the nucleon mass) appears in the numerator of the Feynman amplitude There are 
indeed certain closed three-loop diagrams (see Fig. 3) which contribute to the energy per particle 
already at order O(kj) |T2[. By the same argument one expects that the four-loop Fock-diagramQ 
which includes the twice-iterated one-pion exchange (proportional to M 2 ) will contribute already 
at order O(kj). Luckily, the analogous four-loop Hartree-diagramP] vanishes when taking the spin- 
or Dirac-trace. Because of the enormous complexity of four-loop diagrams we restrict ourselves to 
the computation of all contributions up-to-and-including order O(kj) as they are given by closed 
three-loop diagrams. 

The only new ingredient in performing calculations at finite nucleon density (as compared 
to calculations of scattering processes in the vacuum) is the in-medium nucleon propagator. It 
expresses the fact that the ground-state of the system has changed from an empty vacuum to a 
filled Fermi sea of nucleons. In its relativistic form this in-medium propagator of a nucleon with 

2 Obtaincd by drawing in the right diagram in Fig. 3, a third pion-line intersecting both the other two. 
3 Obtained by drawing in the left diagram in Fig. 3 a third parallel pion-line connecting both equally oriented 
nucleon rings. 
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four-momentum p^ = (po,p) reads 



V + M ){ p 2_M2 + ie - 2 ^^ 2 - M ") *(p<0 0(k f - \p\)) • (3) 

Note that eq.(3) splits additively into the vacuum nucleon propagator and a medium insertion 
(i.e. the piece involving an on-shell delta-function and step-functions). This allows to organize 
the diagrammatic calculation according to the number of medium insertions. Diagrams with no 
medium insertion lead to an unobservable shift of the vacuum energy. Diagrams with exactly 
one medium insertion just renormalize the nucleon mass to its measured value M (because of 
the presence of one single on-shell delta-function). The really interesting many-body effects from 
interactions start thus with diagrams having two or more medium insertions. 

The pion-nucleon interaction vertices relevant in this work are the pseudovector 7riViV-vertex 
and the Tomozawa- Weinberg 7T7riViV-contact vertex of the form, 

77- ialh Ta , J75 _ *U 6 abcT c ■ (4) 

Here, the pion four-momenta q a ^ are out-going ones and f n = 92.4 MeV denotes the weak pion 
decay constant. For the nucleon axial vector coupling constant g& we choose the value qa = 1.3. 
Via the Goldberger-Treiman relation this corresponds to a 7riViV-coupling constant of g^N — 
g^M/ f m = 13.2 which agrees with present empirical determinations of g n N from 7riV-dispersion 
relation analyses |T^] . 



2.1 Kinetic energy 

The first contribution to the energy per particle E(kf) is the kinetic energy of a non-interacting 
relativistic Fermi gas of nucleons. Expanding a/ M 2 + p 2 — M in powers of 1/M and integrating 
over a Fermi sphere of radius kf one gets, 

Ek{kf) = WM~5QM^- (5) 

For the densities p = 2k^/3ir 2 of interest this series converges very rapidly. The next term in this 
series, fc^/48M 5 , is already negligibly small. 



2.2 One-pion exchange 




Fig. 2: The one-pion exchange Fock- diagram. Solid and dashed lines represent nucleons and pions, 
respectively. The combinatoric factor of this diagram is 1/2. 

The one-pion exchange Hartree-diagram is trivially zero (since g M = 0) and the one-pion 
exchange Fock-diagram is shown in Fig. 2. The interesting many-body effect comes from the 



medium insertion at each nucleon propagator. In that case the two loop-integrations convert into 
integrals over Fermi spheres of radius kf. After taking the Dirac-trace of the product of (J>+ M)- 
factors and pseudovector irN- interaction vertices one performs a (non-relativistic) 1/M-expansion 
of the complete integrand. In this form all integrals can be solved analytically and the result for 
the l7r-exchange Fock-diagram reads, 



Ei{k f ) 



3g 2 A ml u 3 1 3u , 3 

, A „ <^ h arctan2w - ( — 
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32m 3 
5) arctan2u 



ln(l + 4w 2 
1 



8m 3 



ln(l + Au 2 



(6) 



Here, we have introduced the abbreviation u = kf/m w . The terms proportional to u 3 and 
u 5 m 2 /M 2 in eq.(6) correspond to the (repulsive) zero-range contact interaction produced by 
any pseudoscalar meson exchange. The 1/M-expansion is again converging rapidly. For example, 
at kf = 270 MeV the l/M 2 -correction (second line in eq.(6)) is only —5.2% of the leading order 
term (first line in eq.(6)). Note also, the Taylor-series expansion of E\(kf) in eq.(6) converges 
only for kf < m n /2. This would correspond to tiny densities of p < 0.0027 fm -3 . 



2.3 Iterated one-pion exchange 




Fig. 3: Iterated one-pion exchange Hartree and Fock- diagrams. The combinatoric factor of these 
diagrams is 1/4- 

Next, we consider the three- loop diagrams shown in Fig. 3. With two medium insertions on 
neighboring nucleon propagators they contribute, via mass and coupling constant renormalization, 
to the one-pion exchange Fock-diagram. These effects are automatically taken care of in eq.(6) 
by using the physical input parameters. Considering two medium insertions on non-neighboring 
nucleon propagators we encounter the planar box graph 0. In a non-relativistic 1/M-expansion 
the first contribution from the planar box graph is the iterated l7r-exchange which is enhanced 
by the large scale factor M. In case of the left diagram in Fig. 3 (Hartree-diagram) the inner loop 
integral becomes equal to the iterated l7r-exchange amplitude in forward direction. Using the 
analytical results given in section 4.3 of ref.0, we can even perform the remaining integral over 
the two Fermi spheres of radius kf. Altogether, we find for the Hartree-diagram in Fig. 3 with 
two medium insertions, 

E 2 (k f ) = [- - 59n + (60 + 32u 2 ) arctan2 M - (— + -) ln(l + An 2 )] . (7) 

Note that this expression does not include the contribution of a linear divergence dl 1 of 
the momentum space loop integral occurring in the Hartree-diagram. In dimensional regu- 
larization such a linear divergence is set to zero, whereas in cut-off regularization it equals 



the momentum space cut-off A. There is an interpretational problem with the result eq.(7) 
when using dimensional regularization. While eq.(7) alone gives sizeable repulsion, for example 
E 2 (kf = 270 MeV) = 50.0 MeV, second order perturbation theory arguments suggest attraction 
(because the intermediate-states lie mostly higher in energy than the initial or final state). The 
source of this apparent contradiction is the linear divergence J °° dl 1. Using cut-off regularization 
together with a large enough value of A one regains the expected attraction. The additional ad- 
ditive terms specific for cut-off regularization are collected in eq.(15). The evaluation of the right 
diagram in Fig. 3 (Fock- diagram) with two medium insertions is very similar. There one encoun- 
ters the iterated l7r-exchange amplitude in backward direction J7|. Putting all pieces together 
one obtains the following result for the Fock-diagram in Fig. 3 with two medium insertions, 



- ., , g\Mmt(u 3 f u , 3x(u - x) 2 (2u + x) 

E 3( k f) = T+ dx ' 



(4tt) 3 /4 t 2 Jo 2u 3 (l + 2x 



2\ 



'l + 8x 2 + 8x 4 ) arctanx — (l + 4x 2 ) arctan2x 



where we have again transferred the linear divergence proportional to the cut-off A to eq.(15). 

Next, we consider the diagrams in Fig. 3 with three medium insertions. In this case one has 
to evaluate an integral over the product of three Fermi spheres of radius kf. Suitable techniques 
allow to perform most of these nine integrations analytically. For the Hartree-diagram in Fig. 3 
with three medium insertions we end up with the following representation, 

Mkf) = 9 £^s £ d ™ 2 f_dy[2uxy + (it 2 - x 2 y 2 )H] { *£±£ - 2 ln(l + s 2 )] , (9) 



H = \ia- — — , s = xy + Ju 2 — x 2 + x 2 y 2 . (10) 
u — xy v 

Similarly, we find for the Fock-diagram in Fig. 3 with three medium insertions the representation, 
9g\Mmj r* \G 2 x 2 f f yz8(y 2 + z 2 - 1) r 2 , . , 2 i 1 

(11) 

with the auxiliary functions, 

I I _)_ (r,, _|_ T )2 

G = U {1 + u* + x 2 ) - — [1 + (u + x) 2 ] [l + (u- x) 2 } In -Z + , (12) 

Ax 1 + (u — xY 



xz 



+ Vu 2 — x 2 + x 2 z 2 . (13) 



For the numerical evaluation of the dydz-double integral in eq.(ll) it is advantageous to first anti- 
symmetrize the integrand both in y and z and then to substitute z = \/y 2 C 2 + 1 — y 2 . This way 
the integration region becomes equal to the unit-square < y, £ < 1. We note that our results, 
eqs.(7,8,9,ll), which comprise all 27r-exchange contributions at order O(kj), agree with those 
of ref. []T2j after suitable simplification. As a further check on the whole formalism we rederived 
the analytical result of Onsager et al. (I7j for the two-photon exchange Fock-diagram of the 
electron-gas, E = [(In 2/3) — (3C(3)/27r 2 )] Rydberg, using the same techniques. 

The diagrams in Fig. 3 with four medium insertions do not have to be calculated explicitly. 
These are purely imaginary and they are cancelled by the imaginary parts coming along with the 
diagrams with two and three medium insertions. Via unitarity, such imaginary parts are related 
to on-shell scattering of nucleons. Energy conservation and Pauli-blo eking in the medium leave 
no phase-space for these processes and thus the energy per particle E(kf) remains a real quantity 
(as it must of course be). 
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2.4 Irreducible two-pion exchange 




Fig. 4- Irreducible two-pion exchange Fock- diagrams. The corresponding Hartree- diagrams vanish 
for isospin symmetric nuclear matter. The combinatoric factors of these diagrams are 1/4, 1, 
1/2 and 1/4, in the order shown. 

Contributions of three-loop 2-7r-exchange diagrams to E(k f ) at order 0(kj) are in accordance 
with the naive chiral power counting rules. The relevant 27r-exchange Fock-diagrams are shown 
in Fig. 4. Considering two medium insertions at the upper and lower nucleon propagators one 
finds that their contribution to the energy per particle can be expressed in terms of the T-matrix 
related to the irreducible 27r-exchange. This T-matrix has been evaluated in section 4.2 of ref.0 
using dimensional regularization. The result of cut-off regularization differs only little by some 
coefficients of purely polynomial terms and the presence of a quadratic divergence A 2 (see also 
appendix A in ref.|J). After performing the remaining integral over two Fermi spheres of radius 
kf one ends up with the following analytical expression for the Fock-diagrams in Fig. 4, 



E 6 (k f ) 



— (43^ + 6g 2 A - 1) + ^(23^ + 2&-1 



\n 2 (u + \/TT 



+ 



XL XL 1 

-(llg 4 A - 10g 2 A - 1) + -(59g A - 50g 2 A - 9) + -(883^ - 90g 2 A - 73) 



+—(1-6^-43^) 



16m 2 
u 



Vl + u 2 \n(u + VI + v 2 ) + —(43^1 + 6g 2 A -l 



32u 



+—(397 + 210£ 2 - 5647^) + -(4 + 5g 2 A + 31g\) + —(119 + 710^ 2 - 349^) 

,5 



+ 



^(15^-6^-1) 



u 



llg\ - 10^1 - 1; 



ln^fl 
2A J 



(14) 



where we have again transferred the term quadratic in the cut-off A to eq.(15). Obviously, terms 
proportional to g A , g A , g A in eq.(14) belong to the first, second, third and fourth diagram in Fig. 4, 
respectively. Note that the last diagram contributes both at order O(kj) and at order O(k^). 
The contribution of the Hartree-diagrams related to irreducible 27r-exchange vanishes in isospin 
symmetric nuclear matter for the following simple reason. After taking spin- and isospin-traces of 
the T-matrix, only the isoscalar central NN-amplitude Vc(0) survives, and this amplitude receives 
no contribution from irreducible chiral 27r-exchange at leading order as shown in section 4.2 of 
ref. J7J. The diagrams in Fig. 4 with three medium insertions contribute only via relativistic 1/M- 
corrections at order 0(k 6 j-). Besides eq.(14), there is no further contribution at order O(kj) from 
three-loop 27r-exchange diagrams. Note also that all contributions to E{kf) remain finite in the 
chiral limit m n = 0. For example, eq.(14) becomes proportional to fc^[ln(&//A) + const]. As far 
as we can see, there are no infra-red singularities generated by the Goldstone bosons (massless 
pions) . 



Finally, we give the complete expression for the power divergences specific to cut-off regular- 
ization, 



A k 



Mkf) = [ ~ Wg\M + (3g 2 A + l)(g A - l)Aj . (15) 

Obviously, the first term proportional to M stems from iterated l7r-exchange. The Hartree and 
Fock-diagrams in Fig. 3 contribute to it in the ratio 4 : 1 as a result of their spin- and isospin 
factors. The term quadratic in the cut-off A, on the other hand, originates from irreducible 
27r-exchange. 

Eq.(15) is equivalent to the contribution of a zero-range NN-contact interaction. The cut-off 
scale A is effectively parametrizing its strength. Note that eq.(15) leads to a sizeable attraction in 
the hundred MeV range for physically reasonable values of the cutoff, 0.5 GeV < A < 1.0 GeV. We 
are thus lead to the conclusion that in isospin symmetric nuclear matter 27r-exchange produces 
"zero-range" attraction and finite-range repulsion. It should also be noted that the dominant (real 
part of) iterated l7r-exchange does not have an equivalent representation as a local coordinate- 
space NN-potential (see section 4.3 in ref.fTj). Such a connection exists only for l7r-exchange and 
(with certain restrictions concerning the zero-range part) for irreducible 27r-exchange |7]]. Note 
that Ek(kf) in eq.(15) encodes the short-distance NN-dynamics necessary for nuclear binding. 
We emphasize that, for reasons of consistency, this piece (which originates from the few divergent 
parts of chiral 27r-exchange) should not by further iterated either with itself or with l7r-exchange, 
in contrast to what has been done with the NN-contact interaction in ref. [12|. 



2.5 Results 

Before presenting results for the nuclear matter equation of state, we would like to exhibit the 
saturation mechanism. For that purpose we truncate the previous calculation of one- and two- 
pion exchange diagrams at order O(kj) and we take the chiral limit, m n = 0, (disregarding for 
the moment the less important A 2 -term in eq.(15)). In that case the complete expression for the 
energy per particle E{kj) turns into the form, eq.(l), with coefficients a and (3 given by 

, , , / 9irN \ 4 A /gnN\ 2 

a 



10 Ii^m-Ut)' (16) 



P = )V 2 + 237 - 24 In 2) - A = 13.55 . (17) 

Note that the parameterfree expression for (3 gives a number quite close to (3 = 12.22 as extracted 
from the realistic equation of the state of ref.[|15[]. Furthermore, by adjusting the cut-off scale 



to A = (0.5. . . 0.6)M (which in fact lies in the physically reasonable range) a will take on its 
required value. Therefore, as long as the effects due to the finite pion mass and the terms at 
order O(kj) do not change this picture completely, realistic nuclear binding is guaranteed by 
chiral pion-nucleon dynamics, together with fine-tuning of the single scale A representing the 
short-distance NN-dynamics. 

Now we include all calculated terms up to order 0(kj) and we set = 135 MeV (the neutral 
pion mass). We fix the (negative) binding energy to the value E = —15.26 MeV obtained in 



extensive and elaborate fits of nuclide masses in ref. |Tj|. The energy per particle E{kf) has 
a minimum with this value of Eq at a density of po — 0.178 fm~ 3 (corresponding to a Fermi 
momentum of k/o = 272.7 MeV = 1.382 fm _1 ) if the cut-off scale is fine-tuned to A = 646.3 MeV. 
Note that a saturation density of po = 0.178 fm -3 is somewhat on the large side. A slight increase 
of gA by 2% would even shift p to the center of the empirical saturation density p — 0.166 fm~ 3 
[O. Let us also investigate how the value E = —15.26 MeV arises in the present calculation. 



Its decomposition into contributions from the kinetic energy and the three classes of diagrams is 
Eq = (23.40 + 18.24 — 68.35 + 11.45) MeV. On the other hand, when ordering in chiral powers 
0(k" f ), = 2, 3, 4, 5) one finds E = (23.75 - 154.54 + 124.61 - 9.08) MeV. The individual entries 
follow closely the pattern prescribed by the parametrization eq.(l) with a balance between large 
third and fourth order terms. If the behavior of the last two entries is representative the chiral 
expansion of E appears to converge. 




Fig. 5: Energy per particle of isospin symmetric nuclear matter derived from chiral one-and two- 
pion exchange (solid line). The value of the cut-off scale is A = 646 MeV. The dashed line 
corresponds to the result of ref. ji^/. 

In Fig. 5 we show the resulting nuclear equation of state for densities up to p = 0.5 fm -3 (i.e. 
about 3po). The nuclear compressibility K related to the curvature of the saturation curve at 
its minimum comes out as K = 255 MeV, in very good agreement with the nowadays accepted 
empirical value K = (250 ± 25) MeV 0, 0] (for the definition of K, see eq.(2)). 

Fig. 6 shows by the solid line the dependence of the saturation point (po, Eq) on the cut-off A 
which has been varied in the range 0.6 GeV < A < 0.7 GeV. The inserted rectangle corresponds 
to the empirical saturation "point" E = (—16 ± l)MeV and kf = (1.35 ± 0.05) fm -1 quoted 
in ref. The variation of the saturation point (po,E ) in Fig. 6 is somewhat reminiscent of 
the familiar Coester-line; of course, the physics behind it is quite different here. It is gratifying 
that our line meets the empirical saturation "point" . The dashed line in Fig. 6, which gives the 
calculated saturation point in the chiral limit, m n = (keeping g^, M, f n fixed), is also interesting. 
It suggests that in QCD explicit chiral symmetry breaking is not a crucial condition for nuclear 
binding. Note however that for fixed A, taking the chiral limit m n = 0, increases the binding 
energy E as well as the equilibrium density p considerably. 
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Fig. 6: The nuclear matter saturation point (p ,E ) at finite pion mass (solid line) and in the 
chiral limit (dashed line) as a function of the cut-off scale A (points in MeV). The inserted 
rectangle corresponds to the empirical saturation "point". 



3 Asymmetry energy 

The properties of isospin symmetric matter can be well reproduced in our framework with just 
one single adjustable cut-off scale A. In order to test the isospin dependence of the underlying 
pion-exchange dynamics which causes saturation, we now investigate the asymmetry energy. In 
isospin asymmetric matter the Fermi seas of protons and neutrons are filled unequally. With the 
help of the projection operators (1 ±r 3 )/2 such an isospin-asymmetric situation is realized by the 
simple substitution, 

0(k f -\p\) - l ^6{k p -\p\)+ l -^6{k n -\p\), (18) 

in the in-medium nucleon propagator eq.(3). Here k p and k n denote the (different) Fermi momenta 
of protons and neutrons. Choosing k PiU = kf(l =F 5) 1 / 3 (with 5 a small parameter) the nucleon 
density p = p p + p n = (k p + kfJ/3-ir 2 = 2k 3 j/3rc 2 stays constant. The expansion of the energy per 
particle of isospin asymmetric nuclear matter, 

E as (k p , k n ) = E(k f ) + 5 2 A(k f ) + ... , (19) 

around the symmetry line (k p = k n or 5 = 0) defines the asymmetry energy A(kf). Note that 
the parameter 5 is equal to (p„ — p p )/(p n + P P ) or (N — Z)/(N + Z). Differences in comparison 
with the diagrammatic calculation in section 2 occur only with respect to isospin factors and the 
radii of the Fermi spheres, k p = kf(l — 5) 1 / 3 or k n = kf(l + 5) 1 / 3 . Following the scheme in section 
2, we summarize the individual contributions to the asymmetry energy A(kf) without going into 
further technical details. 

i) Kinetic energy of a relativistic Fermi gas: 

*<*'> = ™-mfJ- (20) 



1 1 



The next term in this series, k^/16M 5 , is negligibly small. 

ii) l7r-exchange Fock-diagram including the relativistic l/M 2 -correction: 



3 2 

u u 

2~Y + JP 



with the abbreviation u = kf/m n . 

iii) Hartree-diagram in Fig. 3 with two medium insertions: 

g\MmU r 25 7 



u 3 -— - arctan 2u- — - ln(l+4w 2 ) 



(21) 



Mkf) = ft^{(f U + L) ln(1 + ^ ~Y U - 16 ^ arctan2w } • ( 22 ) 



iv) Fock-diagram in Fig. 3 with two medium insertions: 



g\Mm A ( 5 . 



3x 2 - Au 2 



o 6w(l + 2a; 2 ) 



(1 + 8a; 2 + 8a; ) arctan x — (1 + 4a; 2 ) arctan 2a; 



(23) 



v) Hartree-diagram in Fig. 3 with three medium insertions: 



Mkf) = 



g\Mm\ 

(47T/ 7r )%3 

r2s 2 + s A 



X 



x 



1 + s 2 

2 



J dxx 2 J dy 
- 21n(l + s 2 



uxy(26u 2 — 30x 2 y 2 ) 



3(u 2 — x 2 y 2 ) 
Au 2 H s 6 

+ 



+ (3m 2 - 5x 2 y 2 )H 



3(1 + s 2 ) 2 

(5 + s 2 )(3s 2 - 8ss' + 8s' 2 ) + 8s(l + s 2 )(s" - 5s' + 3s) 



2uxy + (u 2 — x 2 y 2 )H 

„4 



3(1 + s 2 ) 3 



(24) 



The auxiliary functions H and s have been defined in eq.(10) and s' and s" denote partial deriva- 
tives, 



ds 



s = u- 



,d 2 



s = u 



(25) 



du ' du 2 
vi) Fock-diagram in Fig. 3 with three medium insertions: 

A 5 (k f ) = j^™l 3 £ dx {^ 3G ™ ~ 2Gu + 3G ° 2 " 8Go1 " 3G) + Ul0 2 i 4 " J ° 1 (3G 10 - 5G i) 



G w + G 



+ X - [ I dyl l dz yz6{y2 + Z2 - 1) 



+ 



2s 3 t 4 (8s' - 3s) 

6" i-i" y J-i^\yz\^y I + z 2 - 1 [(1 + s 2 )(l + t 2 ) 
(3 + s 2 )(8ss' - 3s 2 - 8s' 2 ) + 4s(l + s 2 )(6s' - 3s - 2s") 



s 2 [t 2 -ln(l + t 



2Mn 



(1+S 2 ) 2 



.(26) 



The auxiliary functions G and t have been defined in eqs. (12,13) and we have introduced the 
following double-index notation for partial derivatives of the function G, 



■ ■ d i+j G 

vii) Irreducible 27r-exchange Hartree- and Fock-diagrams in Fig. 4: 

1 



(27) 



Mkf) 



I2u 



'it 

1 - 2g\ - 23g\) + -(l + 2g 2 A - 7g\) \n 2 (u + Vl + u 2 ) 
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+ 



1(23^ + 2^ - 1) + u\\ + 2^ - 5g\) - 



u 



23g 



u 
16 



u 



Vl + u 2 ln(w + VI + u 2 ) 



11 



245 



4 + Y ^ A " "l2~^ 



^(1 + 10^ + ^) 



- + 10g 2 A - 2hg\ - -g\u 2 



2A J 



(28) 



The Hartree-diagrams (not shown in Fig. 4) give a contribution to the asymmetry energy A{kf) 
through the isovector central NN-amplitude Wc(0) (see section 4.2 in ref.0). In terms of 
this amplitude the contribution of the irreducible 27r-exchange Hartree-diagrams reads A{kf) = 
-pW c (0)/2. 

viii) Power divergences specific for cut-off regularization: 



A A (k f ) 



Akj 



3(4tt/ ot 



2Gg 4 A M + 5(3g 2 A + l)(l-g 2 A )A 



(29) 



Again, this expression collects the contributions from iterated l7r-exchange and irreducible 2tc- 
exchange. In comparison to E\(kf) given in eq.(15) the term, eq.(29), has changed sign. The 
isospin-structure of the underlying 27r-exchange interaction shows up in the linear and quadratic 
term in A through relative factors —13/15 and —5/3, respectively. 

As a side remark we note that in the chiral limit, m Tl 
asymmetry energy A(kf) (see eqs. (22,23,24,26)) can be evaluated in closed form with the result 



0, the order 0(kj) contributions to the 



A{k f )\, 



M 



(9Ak f \ 
Utt/J 135 



(447r 2 + 987 + 216hi2) 



(30) 



For the numerical evaluation of the asymmetry energy we use the same parameter input as 
in section 2.5, in particular a cut-off scale of A = 646.3 MeV. The value of the asymmetry energy 
A(kf) at our saturation point kfo = 272.7 MeV comes out as Aq = A{kfo) = 33.8 MeV. This is in 
very good agreement with the empirical value of Aq = 33.2 MeV obtained in ref. [rj| from extensive 
and elaborate fits to nuclide masses. Let us again investigate how the value A = 33.8 MeV arises 
in the present calculation. Its decomposition into contributions from the kinetic energy and the 
three classes of diagrams is A = (12.7 — 5.3 + 40.7 — 14.3) MeV. On the other hand, when ordering 
in chiral powers Oik 



fh 



= 2,3,4,5) one finds A = (13.2 + 129.9 - 127.4 + 18.1) MeV. Here, 
one observes an almost complete cancelation of large third and fourth order terms. Again, if the 
behavior of the last two entries is representative, the chiral expansion of the asymmetry energy 
A should converge. 

In Fig. 7 we show by the solid line the dependence of the asymmetry energy A(kj) on the 
density p = 2A;|/37r 2 . One observes that A(kf) reaches its maximum close to the saturation density 
Pq. We obtain therefore a small value of the (empirically not well determined) slope parameter 
L = kfo[dA(kf)/dkf]o = 14.7MeV (see also ref. ||13[| ). Furthermore, we extract an asymmetry 
compressibility of K as = k 2 Q [d 2 A{kf)/dky\Q — 2L = —359 MeV. Such large and negative values 
of K as are also found in some other calculations (see table 4 in ref.|[L3|). For comparison we 
show in Fig. 7 also the result for the asymmetry energy A(kf) obtained in the Briickner-Hartree- 
Fock calculation of ref.||19||. The downward bending of our asymmetry energy A(kf) at densities 
p > 0.2 fm -3 presumably indicates the limits of validity of the present chiral perturbation theory 
calculation of nuclear matter. 
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Fig. 7: The density dependence of the asymmetry energy A(kf) (solid line). The cut-off' is A = 
646 Me V. The dashed line gives the result of ref. /|7^/. 

The simultaneous correct description of several nuclear matter properties (saturation point, 
compressibility, asymmetry energy) with only one adjustable scale parameter, the cut-off A, is 
indeed highly non-trivial. This strongly indicates that nuclear binding and saturation through 
the combination of chiral one- and two-pion exchange and short- distance dynamics encoded in 
a single high momentum scale A is a realistic scenario (at low densities p < 0.2 fm -3 ). Our 
treatment of the zero-range components of the effective NN-interaction in terms of a cut-off A 
seems to work surprisingly well and it resembles effects of correlation functions in sophisticated 
many-body calculations. 



4 Neutron matter 

The extreme of asymmetric nuclear matter is pure neutron matter. All existing realistic cal- 



culations fl5|, |2(J, g2j agree that pure neutron matter is unbound. Its energy per particle 
monotonically rises with the neutron density. 

In order to arrive at the energy per particle, E n (k n ), of neutron matter in our diagrammatic 
framework it is sufficient to make the substitution 

e{k f -\p\) - l -^e{K-W\), (si) 

in the in-medium nucleon propagator eq.(3). Here k n denotes the Fermi momentum of the neu- 
trons related to the neutron density by p n = /c 3 /37r 2 . As a consequence of the substitution eq.(31) 
only the isospin factors of individual diagrams change and all Fermi spheres have the radius k n . 
We follow the scheme in section 3 and enumerate the individual contributions to the energy per 
particle, E n (k n ), of neutron matter. 

i) Kinetic energy of a relativistic Fermi gas: 

E nk (k n ) = ^-^^. (32) 
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ii) l7r-exchange Fock-diagram including the relativistic l/M 2 -correction: 



E nA (k n ) 



g A m% f m 



(47T/ ff )2 I 3 
,2 r 40 



5// 



3u 
T 



3 1 

^ + 7T- - + arctan 2m - ( — + — -) ln(l + 4m 2 

Hii A/hp I 



ml 



40M 2 



— m 3 - 8u 5 + 9u + — - (12m 2 + 5) arctan 2u — ln(l + 4m 2 ] 

3 2m 8u 3 



(33) 



We emphasize that in this section the meaning of u changes to u — kn/m^. 
iii) Hartree-diagram in Fig. 3 with two medium insertions: 



- g\Mm\ f 9 r 
E n:2 {k n ) = - < 59m 



10(8tt) 3 /4 12m 



(60 + 32m 2 ) arctan 2m - f + — ) ln(l + 4m 2 ) } . (34) 

v 8m j 2m 7 J 



iv) Fock-diagram in Fig. 3 with two medium insertions: 



E n $(k r , 



g\Mm% f M 3 _f u ^ x(u - x) 2 (2m + x) 



(4tt) 3 /4 [ 6 Jo 2u 3 (l + 2x 2 ) 
v) Hartree-diagram in Fig. 3 with three medium insertions: 



l+4x ) arctan 2x— (l+8a; +8x )arctanx 



(35) 



Sn,^n = f. A . u 3 / cfex 2 / 2mxj/ + m 2 - x 2 y 2 )H 

(47T /tt) M J J-l L J 12(1 + 5^) 



ln(l 



(36) 



vi) Fock-diagram in Fig. 3 with three medium insertions: 



(4tt/,)% 3 A 



4 J-i V-i Vz/ 2 + 2 2 — 1 



s 2 -ln(l+s 2 ) t 2 -ln(l+t 2 

(37) 



vii) Irreducible 27r-exchange Hartree- and Fock-diagrams in Fig. 4: 

,5 ,r 1 I 



E n fi(k Ti 



(47T/.) 4 
+ 



, 2//3 v834 + 6<? 2 - 1) + ^(47^ + 2^ 2 - i; 



ln 2 (u + Vl +m 2 ) 



-(34 - 50<& - 9) - —(13^ + 10^ 2 + 1) + j2q(1691^ - 90g 2 A - 73) 



1 



16u 2 
u 



[l-6g 2 A - 83g\) 



VTT^Hu + VTT^) + ^ r (83g A + Qg 2 A - 1) 

32m 



u 



480 

u 3 



(397 + 210g% - 11159^) + (119 + 710^ + 107# 



1800 



(3* 



+ ^(9 + 55<? 2 -108^)4- 



u 



■'(- + 2g\ - 5g A ) - ^(13g A + 10g 2 A + 1) 



In^l 
2A J 



This expression includes the contributions from Hartree diagrams (not shown in Fig. 4) via the 
isovector central NN-amplitude Wc{0) (see section 4.2 in ref.Oj). In terms of this amplitude the 
contribution of the irreducible 27r-exchange Hartree-diagrams is, E n {k n ) = — p n Wc(0)/2. 



viii) Power divergences specific for cut-off regularization: 



Akl 



3(47r/ ff ; 



2g\M+(3g 2 A +l){g 2 A -l)A 



(39) 



This 



Note that compared to eq.(15) the attraction in neutron matter gets strongly reduced, 
welcome feature has its origin in the isospin dependence of the 27r-exchange. 

In the chiral limit, m n = 0, and truncated at order O(kj) the equation of state of neutron 
matter E n (k n ) turns also into the form eq.(l) with the coefficient (3 n = (g^Ni^n) 4 (27 — Air 2 — 
241n2)/70 — (3/56) = —0.562. The other coefficient a n can be easily read off from eqs.(33,39). 




Fig. 8: The equation of state of pure neutron matter E n (k n ). The dashed line corresponds to 
the many-body calculation of ref. [T3j]. The dashed-dotted line represents eq.(40) with a neutron- 
neutron scattering length of a nn = 18.7 fm. The full line gives the result of chiral one- and two-pion 
exchange with a cutoff of A = 646 Me V. 

If we use the same parameter input as in section 2.5, the resulting equation of state of neutron 
matter, E n (k n ), comes out as the solid line in Fig. 8. The convex shape of the curve is generic 
and does not change much with the cut-off A. It demonstrates the limitations of our present 
perturbative approach. The dashed line in Fig. 8 corresponds to the many-body calculation of 
the Urbana group jl5fl . This curve should be considered as a representative of the host of existing 



neutron matter calculations |p0| , [2T| , pqj which scatter around it. One observes a rough agreement 
up to neutron densities of about p n = 0.25 fm~ 3 . At higher densities our neutron equation of state 
starts to become unrealistic because of its downward bending (/3 n < 0). This feature is inherited 
from the asymmetry energy A(kf) (see Fig. 7). Note however that our prediction for the neutron 
matter equation of state is still much better than the one of the cro;-mean field model |23|] (for 
low neutron densities p n < 0.25 fm -3 ). 

The complete resummation of in-medium multi-loop diagrams for a system with an unnaturally 



large scattering length has been achieved in the limit of large space-time dimensions D in ref. [11 
Neglecting the l/D-corrections (approximately less than 20%) the result of Steele [|TTJ applied to 
neutron matter reads, 

k^ f 3 a nil k n 



En(k n ) = ^ 



10 37r + 6a nn k n 



(40) 



with a nn the neutron-neutron scattering length. A recent measurement P3 has found the precise 
value a nn = (18.7 ± 0.6) fm. (We use that sign- convent ion in which a positive scattering length 
corresponds to attraction.) The density dependence of eq.(40) is shown by the dashed-dotted 
line in Fig. 8. The agreement with the many-body calculation of ref. |15[ is somewhat better.^ 
Interestingly, the dashed-dotted line does almost not change if the nn-scattering length is sent to 
infinity, a nn = oo, where E n {k n ) = 2k 2 l /lhM. 

Concerning our neutron matter equation of state (solid line in Fig. 8) one should not forget 
that it has no free parameters (assuming that the high momentum scale A is the same as in isospin 
symmetric nuclear matter). The fair agreement with the other two curves (for p n < 0.25 fm -3 ) 
is therefore quite surprising. The mere fact the neutron matter is predicted to be unbound is 
already non-trivial. 



5 Summary and Outlook 

The present work can be summarized as follows: We have used chiral perturbation theory to 
calculate the nuclear matter equation of state systematically up to three-loops. The contributions 
to the energy per particle E(kf) from one- and two-pion exchange diagrams are ordered in powers 
of the Fermi momentum kf (modulo functions of kf/m^). We have evaluated all contributions 
up-to-and-including order O(k^). 

A momentum cut-off scale A has been used to regularize the few divergent parts associated with 
chiral 27r-exchange. The terms linear and quadratic in A effectively parametrize an equivalent zero- 
range NN-contact interaction which is strongly attractive in isospin symmetric nuclear matter. 
Its isospin dependence is revealed (and tested) in the asymmetry energy A(kf). 

The saturation mechanism behind the chiral 27r-exchange is easy to understand (in the chiral 
limit m-x = 0). An attractive a k^/M 2 term plus a repulsive (3 kj/M 3 term, with proper coefficients 
a and (3, lead automatically to a realistic nuclear matter equation of state. 

Without inclusion of any further short-range terms, and adjusting the cut-off scale to the physi- 
cally sensible value A ~ 0.65 GeV, the empirical saturation point (po — 0.17fm~ 3 , E ~ —15 MeV) 
and the nuclear compressibility K ~ 250 MeV are well reproduced. Decomposing the binding en- 
ergy Eq into chiral orders one recovers the pattern of the realistic (a, /3)-parametrization. We 
note in passing that the present calculation gives a single particle potential for a nucleon at rest 
in equilibrium nuclear matter which amount to U(0, kfo) ~ —53 MeV ||25|| . 

In the same framework and using the same parameters the density dependent asymmetry 
energy A(kf) has been calculated. The prediction for the asymmetry energy at the saturation 
point Ao ~ 34 MeV is in very good agreement with the empirical value. This is achieved again 
with just one cut-off A instead of introducing additional contact terms: apparently, one single cut- 
off scale is enough to represent the relevant short-distance dynamics. The downward bending of 
A(kf) at p > 0.2 fm -3 presumably indicates the limits of validity of the present chiral perturbation 
theory calculation of nuclear matter. 

The equation of state of pure neutron matter as predicted in this framework is also in rough 
agreement with sophisticated many-body calculations and a resummation result of effective field 
theory for neutron densities p n less than 0.25 fm -3 . The mere fact that neutron matter comes out 
unbound is non-trivial. 

The present approach to the nuclear matter problem is quite different from most other com- 
monly used ones. We do neither start from a (so-called) realistic NN-potential which fits the 
deuteron properties and the NN-phase shifts, nor do we employ relativistic mean field methods. 

4 We have assumed here that eq.(40) can be extrapolated to Fermi momenta as large as k n — 385 MeV. 
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Nevertheless, we find realistic nuclear binding, including its isospin dependence, already at the 
level of the three-loop approximation of in-medium chiral perturbation theory. 

Questions remain, of course, whether the good results obtained so far survive in higher orders. 
Relativistic 1/M-corrections seem to present no problem since they are in general small (as 
long as the densities are not too high). Two-pion exchange with virtual A(1232)-resonance 
excitation produces essentially all the medium-range scalar-isoscalar attraction which is needed 
for the description of the peripheral NN-scattering ||. In the nuclear matter calculation the 
short-range contributions from 27r-exchange with A(1232)-excitation as well as additional Pauli- 
blocking effects (related to diagrams with three medium insertions) come also into play. At 
the present stage, such effects are possibly hidden in the cut-off scale A; at a higher level of 
"resolution", they can and have to be studied explicitly in the framework of chiral effective field 
theory. Work along this line is in progress. 
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